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Abstract 
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1 Introduction 



The Skyrme model [8] is a nonlinear theory of pions in which baryons are described by 
soliton solutions. With only pion degrees of freedom the Skyrme model requires the inclusion 
of the Skyrme term, which is quartic in derivatives, in order to yield stable soliton solutions. 
Many years ago it was realized that including u mesons, in addition to pions, produces a 
stable soliton solution without the need for a Skyrme term PQ. Unfortunately, numerical 
computations of multi-solitons in the u meson model is a formidable challenge, and to date 
no multi-soliton solutions have been computed. This is in contrast to the Skyrme model, 
where numerical solutions have been obtained up to reasonably large soliton numbers [3]. 
However, by using analytical approximations, recent progress has been made on the study of 
multi-solitons in the u meson model [H] , and the results suggest a surprising similarity with 
solitons in the Skyrme model. 

The baby Skyrme model [7] is a (2+l)-dimensional analogue of the Skyrme model. It is a 
modified 0(3) a- model that includes a Skyrme term. In this paper we study the soliton solu- 
tions of an alternative to the baby Skyrme model, in which the Skyrme term is removed and 
stabilization is achieved by including a massive vector meson. A comparison is made between 
the multi-solitons in the two models, as this is a more tractable lower-dimensional analogue 
of the (3+l)-dimensional Skyrmion situation discussed above. We find a remarkable simi- 
larity between the multi-solitons in the two models, at both the qualitative and quantitative 
level, and we explain this behaviour by using a simple derivative expansion. These results 
provide a further justification for the approximate techniques used in (3+l)-dimensions [9], 
that motivated the present work. 



2 The baby Skyrme model and vector mesons 

The Lagrangian density of the baby Skyrme model in (2+l)-dimensions is given by [7] 

C BS = l -8^ ■ d»cj> - m 2 (l - 3 ) - ^{d^ x d v <f>) ■ x d», (2.1) 

where <f> = (0i; 02, 03) is a three-component unit vector 0-0 = 1, and m is the mass of the 
0i and 02 fields, which are the analogues of the pions fields in the Skyrme model. The final 
term in (12. ip is the Skyrme term, with k a positive constant. The size of the baby Skyrmion 
is determined by the combination Jk/tti. 

Finite energy requires the boundary condition <fi — > e 3 = (0,0,1) as |x| — > oo, which 
implies that the spatial plane is compactified to the two-sphere. Therefore, at fixed time, 
<p : S 2 i— > S 2 , with an associated integer winding number, B e Z = -k 2 {S 2 ), which is the 
soliton number. 

There is a conserved topological current 

B/« = _I e ^.(9 a X 9^). (2.2) 
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and the soliton number is the integral over space of the topological charge density B°, that 
is, B = J B° d 2 x. 

Baby Skyrmion solutions for 1 < B < 6 were computed numerically in [TJ, using the 
parameter values m = 1/y/lO and k — 1. Briefly, the solitons are all bound states and are 
axially symmetric for B < 2, but have only discrete symmetries for B > 2. We shall discuss 
these solutions in more detail in the following Section. 

Our alternative to the baby Skyrme model is given by the Lagrangian density 

£ VM = ^<M M 0-m 2 (l-03)-|(d^-3^ ( 2 -3) 

where the Skyrme term has been removed and a vector field uj^ with mass M has been added. 
This is the analogue of the u meson in the (3+l)-dimensional theory, hence the notation. 
The coupling of the vector field to the topological current, with positive coupling constant 
g, also mirrors the higher dimensional theory [I]. 

In the remainder of this paper we are concerned with static solutions of the vector meson 
theory ( 12.31) . For static fields the spatial components of the topological current vanish B l = 
and therefore u>i = 0, since the topological current B^ provides the source for uj^. 

Only static fields are considered from now on, so for notational convenience we write 
uj = ujq. With uji = 0, the static energy derived from (12. 3p is given by 

e vm = I Ag i(f) . di(j) + m 2 (1 _ 03) _ \ d . ud .u _ I M V + ^-aje ij( j> ■ (9,0 x d^)) d 2 x. (2.4) 
The static field equations that follow from the variation of (12 .4p are 

didi<p + m 2 e 3 + -jj-eijdjU (f) x 3$ + (<9j0 • - m 2 3 )0 = 0, (2.5) 

47T 

and 

OAoj - M 2 u = - J^;<£ • {di<t> x d j( f>). (2.6) 

For fields that satisfy equation (12. 6p then multiplication of this equation by u and its substi- 
tuion into (12.41) . together with an integration by parts, allows the energy (12.41) to be rewritten 
in the form 

E VM = I . d%4) + m 2 (1 _ + 1, ^ . ((9 . x Q d 2^ (2 J) 

J x 2 l07T ' 

which will be convenient later. 

One way to see how the vector meson theory evades Derrick's theorem [6], on the non- 
existence of solitons, is to formally solve ( 12.6f) for uj in terms of its Green's function and 
substitute this back into the energy expression (12. 7p . This formulation presents the model in 
terms of a non-local interaction, and is similar to planar soliton models used in the study of 
quantum Hall ferromagnets [10]. However, a more informative approach is to approximate 
the solution of (12.61) by applying a derivative expansion. In fact, only the leading order term 
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is required, which simply corresponds to neglecting the Laplacian term in (12. 6p so that the 
approximate solution is simply 

Substituting this approximation into the energy (I2.7p gives 

£VM „ | . ^ + m 2 (1 _ + ^^\d l( f> X «9 2 0| 2 ) rf 2 X, (2.9) 

which is precisely the static energy of the baby Skyrme model (12. ip upon identification of 
the parameter k as 

K = -i—. (2.10) 
4vrM v ; 

This analysis suggests that soliton solutions should be similar in the vector meson and baby 

Skyrme models. In the following Section we shall confirm this expectation by presenting the 

results of numerical computations of solitons. 

Note that the approximation used in the above analysis becomes more accurate as M 

increases, as the mass term in (I2.6P is increasingly dominant over the neglected Laplacian 

term. In particular, the analysis is not valid in the massless case M = 0. In this case equation 

( 12. 6p reveals that uj does not scale with a rescaling of the spatial coordinates, therefore the 

interaction energy in (12. 7p has the same scale invariance as the cr-model energy. A rigorous 

mathematical analysis of the existence and uniqueness properties of theories of this type 

defined on a torus can be found in [5]. In the massless limit M = 0, lump-like solutions 

with an arbitrary scale can only exist if m = 0, so that the total energy is scale invariant. 

It might be interesting to investigate the lump solutions of such a doubly massless model, 

M = m = 0, but we shall not pursue this here. 



3 Soliton solutions 

As in the baby Skyrme model, we expect that for B = 1 and B = 2 the minimal en- 
ergy solitons are axially symmetric. The axially symmetric ansatz has the form u(r), with 
boundary conditions o/(0) = 0, u(oo) = 0, and 

4> = (sin / cos B9, sin /sin B9, cos /), (3.1) 

where the profile function f(r) has boundary conditions /(0) = it and /(oo) = 0. 
With this ansatz the static energy (12. 4p becomes 

p vm o (°°( l f'2 . B 2 sin 2 f 2 1 n 1 2 2 gB ,sin/\ 

E = 2n (n/ +^ 5— + m (l-cos/) --u --M u +— uf )rdr, 3.2 

jo 2 2 r z 2 2 4tt r ' 

and the static field equations (12. 5p and (12. 6p reduce to the two ordinary differential equations 

f + ~f ~ sm ( 2 /) ~ m sm / + -ir u = °» ( 3 - 3 ) 

r 2r z Air r 
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J' + -LU' 

r 



gB sinf 



0. 



(3.4) 



Arc " r 

In order to compare the vector meson model with the baby Skyrme model we need to choose 
some parameter values. For the baby Skyrme model we fix the parameters to those used in 
the original investigations |7J and given earlier as m = l/ylO and k — 1. Using the same 
numerical codes that we apply below to the vector meson model, we have recalculated the 
energies of solitons in the baby Skyrme model. For B = 1 and B = 2 we compute the 
baby Skyrmion energies to be Ef s = 19.66 and Ef s = 36.90, which agree with the values 
presented in [7J to a good accuracy. 

Next we turn to choosing the parameters in the vector meson model. Motivated by the 
higher-dimensional theory, we choose M so that the ratio m/M is of a similar order to the 
ratio of the pion to u meson mass. The value M = 3/2 is reasonable from this point of 
view and we take this from now on. Given the earlier comments, the expectation is that the 
results will not be too sensitive to the value of M, providing it is sufficiently large. 

The only remaining parameter is g. Given the values of the other parameters, the formula 
(I2.10p suggests g = AttMk = 6tt ~ 18.85. However, rather than using this value, which relied 
upon the approximation (12.81) . we fix g by requiring that the energy of the B = 1 soliton is 
the same as in the baby Skyrme model. This results in the slightly larger value g = 20.83, 
which we shall use from now on. 

Solving the ordinary differential equations (13.31) and (13.41) using a heat flow method yields 
the energies E^ M = 19.66 and E™ = 37.32 < 2E™. By construction E™ = Ef s , but it is 
a non-trivial result that E^ M and Ef s are very close: Table [T] lists the energies in the two 
models for ease of comparison. 



B 






1 


19.66 


19.66 


2 


36.90 


37.32 


3 


55.58 


56.19 


4 


73.61 


74.48 



Table 1: Energies of the minimal energy solitons with 1 < B < 4. Here E is the energy 
in the baby Skyrme model and E VM is the energy in the vector meson model. 



Figured] presents the associated profile functions for B = 1 and B = 2, with the two solid 
curves being those of the vector meson model and the two dashed curves those of the baby 
Skyrme model. This demonstrates the remarkable similarity between the soliton solutions 
of the two models, and goes beyond the agreement of the energies. 

As in the baby Skyrme model [7], it turns out that in the vector meson model the minimal 
energy solitons for B > 2 are not axially symmetric. In fact, by computing the energies of 
axially symmetric solutions, it is already clear that there is no axially symmetric bound state 
with B = 4 since it has an energy of 77.63 > 2E^ M . 

In order to study solitons with B > 2, numerical solutions of the full two-dimensional 
static field equations (12.51) and (I2.6P must be computed. The numerical algorithm uses 
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Figure 1: Comparison of the profile functions for B = 1 and B = 2 solitons in the vector 
meson model (solid curves) and the baby Skyrme model (dashed curves). The lower curves 
correspond to B — 1 and the upper curves to B — 2. 
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Figure 2: The axially symmetric field u(r) for solitons in the vector meson model with 
B = 1 (solid curve) and B = 2 (dashed curve). 
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a heat flow method applied to these equations, with spatial derivatives approximated by 
fourth-order accurate finite differences with a lattice spacing Ax = 0.2 and a grid containing 
200 x 200 lattice points. 

For comparison, we also compute the soliton solutions of the baby Skyrme model using 
the same code. For B = 1 and B = 2 the energies calculated using the two-dimensional 
code agree with those of the axially symmetric computations to the accuracy presented in 
Table[TJ Note that our energy values for B > 2 differ slightly from those presented in [7]. 
Comparing with the axially symmetric computations for B = 1 and B = 2, reveals that 
the two-dimensional results of [7\ under-estimate the energy by approximately 1%, whereas 
our values agree to a greater accuracy. This improvement is probably a result of our use of 
fourth-order accurate finite differences, since we also find a similar under-estimate if we use 
only second-order accurate finite differences. 

X A. 

Figure 3: Plots of the topological charge density for soliton solutions with 1 < B < 4. The 
left-hand-side images are for the baby Skyrme model and the right-hand-side images are for 
the vector meson model. It is clear that there is a remarkable similarity between the solitons 
in the two models. 

In the last two rows of Table [I] we list the energies for the B = 3 and B = 4 solitons 
in the two models, computed using the two-dimensional code. The energies in this table 
confirm that all the solitons are bound states. Even for the non-axial solitons, the energies 
in the two models are very similar. In Figure [3] we present plots of the topological charge 
density B° for 1 < B < 4 in both models. The images on the left-hand-side correspond to 
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the baby Skyrme model and those on the right-hand-side to the vector meson model. These 
plots confirm the amazing similarity between the solutions of the two theories. 

A modification of the baby Skyrme model was introduced in [IT], in which the second 
term in (12. ip is replaced by the symmetric mass term 

C sym = - l -m 2 {l-<Pl). (3.5) 

As this term is symmetric under the replacement 03 i— > —03, then the associated contribution 
to the energy density vanishes at the centre of the soliton, where 03 = —1. This contrasts with 
the original asymmetric mass term in (12.11) . which gives an associated maximal contribution 
to the energy density when 03 = —1. Axially symmetric baby Skyrmions with B > 1 have 
an energy density that is maximal on a circle, and may be viewed as circular domain walls 
separating the vacuum 3 = 1 outside the wall from the vacuum 3 = — 1 inside the wall. 
For B > 2 such configurations are not favourable with the asymmetric mass term, because 
this results in a contribution to the energy that is proportional to the area enclosed within 
the wall. However, for the symmetric mass term (13. 5p this contribution is removed and it 
turns out that the axially symmetric soliton has minimal energy for all values of B [IT] . A 
similar phenomenon occurs in the (3+l)-dimensional Skyrme model, where Skyrmions form 
hollow polyhedral shells without a mass term [3], but an asymmetric mass term produces 
more compact configurations that are no longer shell-like [U [2] . 

To further examine the similarity between the Skyrme term and vector meson stabiliza- 
tion, we compute axially symmetric solitons in the baby Skyrme and vector meson models 
with the asymmetric mass term replaced by the symmetric mass term (I3.5P in both theories. 
The corresponding energies for solitons with 1 < B < 4 are presented for both theories 
in Table El Note that even though we retain the value g = 20.83, rather than fixing g by 
equating the single soliton energies in the two theories, all the soliton energies are again very 
close. In particular, the vector meson model has axially symmetric bound states, sharing 
the same properties as the baby Skyrme model with the symmetric mass term. 



B 


E 




1 


18.18 


18.30 


2 


32.91 


33.28 


3 


48.29 


48.82 


4 


63.89 


64.58 



Table 2: Energies of the (axially symmetric) minimal energy solitons with 1 < B < 4, 
using the symmetric mass term. E BS is the energy in the baby Skyrme model and E is 
the energy in the vector meson model. 



4 Conclusion 

We have presented the results of numerical computations of multi-solitons in an alternative 
to the baby Skyrme model, in which the Skyrme term is removed and stabilization is provided 
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by the coupling to a massive vector meson field. These results reveal a remarkable similarity 
between the solitons of the two theories, which can be understood using a simple derivative 
expansion. This is evidence for a similar correspondence in the (3+l)-dimensional Skyrme 
model, which has been studied recently using approximate methods, but has so far not 
received the detailed numerical investigation that we have been able to perform here in this 
lower-dimensional analogue. 

There are several possible avenues for further investigation of the vector meson theory 
introduced in this paper. The vector meson field introduced here is the lower-dimensional 
analogue of the u meson, so a natural extension is to include the lower-dimensional analogue 
of the p meson. Also, all studies in the present paper have been concerned with static 
solitons, and it would be interesting to study soliton dynamics in this theory and make a 
comparison with dynamics in the baby Skyrme model. 
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